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James York, in a major extension of Andre Lichnerowicz's work, showed how to construct so- 
lutions to the constraint equations of general relativity. The York method consists of choosing a 
3-metric on a given manifold; a divergence-free, tracefree (TT) symmetric 2-tensor wrt this metric; 
and a single number, the trace of the extrinsic curvature. One then obtains a quasi-linear elliptic 
equation for a scalar function, the Lichnerowicz- York (L-Y) equation. The solution of this equation 
is used as a conformal factor to transform the data into a set that satisfies the constraints. If the 
manifold is compact and without boundary, one quantity that emerges is the volume of the physical 
space. This article reinterprets the L-Y equation as an eigenvalue equation so as to get a set of data 
with a preset physical volume. One chooses the conformal metric, the TT tensor, and the physical 
volume, while regarding the trace of the extrinsic curvature as a free parameter. The resulting equa- 
tion has extremely nice uniqueness and existence properties. A even more radical approach would 
be to fix the base (conformal) metric, the physical volume, and the trace. One also selects a TT 
tensor, but one is free to multiply it by a constant (unspecified). One then solves the L-Y equation 
as an eigenvalue equation for this constant. A third choice would be to fix the TT tensor and and 
multiply the base metric by a constant. Aech of these three formulations has good uniqueness and 
existence properties. 



I. INTRODUCTION 



When general relativity is considered as a dynamical system, the Einstein equations, just like the Maxwell equations, 
split into constraint and evolution equations. In this article I focus entirely on the constraints. The initial data consists 
of a pair (gij,K v ), defined on a given 3-manifold, where g^ is a Riemannian 3-metric and K % i is a symmetric tensor. 
K lJ is the extrinsic curvature of the 3-manifold when regarded as an embedded hypersurface of a pseudo-Riemannian 
4-manifold that satisfies the Einstein equations. [1] 

The initial data cannot be specified freely: they must satisfy both a scalar and vector equation, known respectively 
as the Hamiltonian and momentum constraints, i.e., 

R-K i:j K i: > +K 2 = 0, (1) 



V i (K i t-g i *K) = 0, (2) 

where R is the scalar curvature and K = gijK 1 ^ is the trace of the extrinsic curvature. This trace K will play a key 
role in the rest of this article. 

I (implicitly) define the extinsic curvature by writing the relation between the extrinsic curvature and the time 
derivative of the metric as 

d t gij = 2NKij + ViNj + , (3) 

where N is the lapse and N l is the shift. This means that my sign convention for agrees with Wald [3] and is the 
opposite of Misner, Thorne, Wheeler [2]. In particular 

^=d t y/g = NK + V i N i . (4) 
V9 

The first significant step in understanding how to solve the constraints was achieved by Andre Lichnerowicz in 1944 
[4]. He realised that if the extrinsic curvature was tracefree, then the momentum constraint, Eq.(2), required that the 
extrinsic curvature be both divergence- free and tracefree (described as transverse tracefree, hence TT). TT tensors 
are conformally covariant. Given a TT tensor wrt a given metric, on making the following transformations 

—> 4>*9H, Kf T = 4>- w K l 4 Tl (5) 
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then K^ T is TT wrt the transformed metric 7^7. Thus it continues to satisfy the tracefree momentum constraint. 
Lichnerowicz had the idea of choosing a base metric and a TT tensor on it, and making a conformal transformation 
so that the transformed metric would satisfy the Hamiltonian constraint, i.e., 



R-K? T K^ = 0. (6) 

This reduces to the well-known Lichnerowicz equation for the conformal factor (j> 

8W 2 4> -R<P + K l i> T K£ T (j>- 7 = 0. (7) 

Of course, since the solution 4> is to be used as a conformal factor, one is only interested in solutions which are 
everywhere positive. 

It turns out that for either the asymptotically flat case or for the compact, without boundary, case, the Lichnerowicz 
equation can only be solved for a restricted class of data. The choice of TT data is essentially irrelevant, but the choice 
of the base metric is important. The 'tracefree' Hamiltonian constraint Eq.(6) requires the solution scalar curvature 
to be everywhere nonnegative. There exists a conformal invariant (the Yamabe constant) [5] 

y = wM?, (8) 
[fe e dv]i 

where the infimum is taken over smooth functions in the compact case and over smooth functions of compact support 
in the asymptotically flat case. The minimizing function solves the Yamabe equation, and, when used as a conformal 
factor, transforms the manifold into a compact one with constant scalar curvature. On such a manifold, the minimizing 
function can be chosen equal to 1, and then we get 

Y =XM=R vl (9) 

A given metric can be conformally transformed into a metric with positive scalar curvature if and only if the Yamabe 
constant is positive. As a result, the Lichnerowicz equation can be solved if and only if the metric belongs to the 
positive Yamabe class. 

In 1970, James York [6] extended the Lichnerowicz approach significantly when he realised that one could add a 
constant trace term, K, to the extrinsic curvature without disturbing the conformal covariance of the momentum 
constraint, i.e., 

07,,,. m = <f>- w K? T + \k^. (10) 

Note that the constant trace term is a conformal invariant! The equation for the conformal factor becomes the 
Lichnerowicz - York equation 

8V 2 - R4> + K l 4 T Kj^- 7 - \k 2 ^ = 0. (11) 

The Hamiltonian constraint no longer places any restriction on the sign of the scalar curvature when the trace of 
the extrinsic curvature is nonzero. Therefore the Lichnerowicz - York equation can almost invariably be solved with 
any choice of base metric, any choice of TT tensor, and any choice of K. It works both in the compact and the 
asymptotically flat cases. The initial data that is constructed has constant trace of the extrinsic curvature and is thus 
called a constant mean curvature (CMC) solution to the constraints. 

The standard method is quite straightforward: one chooses the base metric, a TT tensor, and a value of K ; one then 
solves for the conformal factor and constructs the physical metric; given the physical metric, one can find the physical 
volume. The aim of this article, however, is to show that the 'standard' interpretation of the Lichnerowicz - York 
equation, at least in the compact case, is not the only reasonable one. It turns out that one can trade off knowledge 
of K for a specification of the physical volume (V p = J <fr 6 dv). This means reading Eq.(ll) as an eigenvalue equation 
for K by putting a normalization condition on the conformal factor so that the physical volume has a prespecificd 
value. 

There are even more radical interpretations. If one is given a tensor which is TT wrt some metric, then one can 
multiply it by any constant and it remains TT. Now solve the Lichnerowicz - York equation, specifying the base 
metric, the TT tensor defined up to a multiplicative constant, the value of K, and the physical volume. The freedom 
to choose the factor multiplying the TT tensor means that this is also a wellposed system. 

Finally, one could fix the TT tensor, but allow the base metric be defined up to a multiplicative constant. Again, 
both K and V can be fixed and the constant regarded as an eigenvalue. These three options will be discussed further. 
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II. TRADING OFF K FOR V P 



As part of his analysis of the constraints, James York [6] suggested that conformal superspace - the space of 
Riemannian 3-metrics, modulo diffcomorphisms, and conformal transformations [8] - should be regarded as the natural 
configuration space of GR. Recently, a number of us [9] have found an action on conformal superspace which generates 
GR in the CMC gauge. A major difference here is that we need to specify the physical volume as part of our initial 
data, while the value of K emerges as part of the solution. An obvious question is whether this switch loses us the 
extremely nice existence and uniqueness properties of the original form of the Lichnerowicz - York equation. The 
answer, as revealed here, is that we lose essentially nothing. 

As can be seen from Eq.(4), and as previously pointed out by York [6], K and the volume can be regarded as 
canonically conjugate variables; therefore switching from K to V p should be viewed as a Legandre transformation, 
and so should not cause any major disruption of the system. In [9] we write down an action, so we are working in 
the Lagrangian rather that the Hamiltonian framework. Therefore it seems more natural to use the volume (a metric 
quantity) rather than specifying K, which is part of the conjugate momentum. Nevertheless, these sorts of arguments 
can only be regarded as poor substitutes for showing that the equation itself is well behaved after the switch from K 
to V p . 

What we need to do is to fix the (conformal) metric and the TT tensor, and to allow K 2 to span the entire range 
from to oo, solving the Lichnerowicz - York equation for each value of K 2 . We should then investigate what happens 
to the physical volume as a function of K 2 . Note that the sign of K does not matter in the Lichnerowicz - York 
equation: only K 2 appears. The sign plays a role only in the evolution, when we determine whether the universe 
either expands or contracts to the future. Let us start with a perturbation calculation. 

We begin with a solution to the constraints, i.e., we are given a metric, a TT tensor, and a constant K, which 
satisfy 



We perturb this data by changing only K 2 , by an amount SK 2 , and solve the perturbed Lichnerowicz - York equation 
(remember we are perturbing around </>=!) 



This is an inhomogeneous linear elliptic equation for 5<j). We can simplify it by eliminating R using Eq.(12) to get 



This is a particularly nice elliptic equation, because the coefficient of the linear term is everywhere negative, which, 
in turn, guarantees that a unique solution to the inhomogeneous equation exists. This should not come as a surprise, 
because we know that the nonlinear equation always has a regular solution. What is interesting, however, is that we 
can use the maximum principle to control the sign of 5<j). Let us assume that SK 2 is positive. If 8<p were positive 
in any region, we would have a point where 5(f) has a positive maximum. However, at a positive maximum, \7 2 8<p 
is nonpositive, and the linear term is negative, which is incompatible with a positive right hand side. Therefore, we 
know that 5<p < 0. This means that as K 2 increases (j> decreases at each point, and the physical volume monotonically 
decreases. The mapping between V p and K 2 is one to one. 

We have a fixed base metric, a fixed TT tensor, and a K 2 which can become either unboundedly large or go to zero. 
The next question to ask is what happens when K 2 becomes large. It can be shown that (j) shrinks monotonically to 
everywhere, and the physical volume vanishes. The easiest way to see this is to consider the alternative. For each 
choice of K 2 we solve for </>, which we know decreases. Let us assume that this 4> does not go to zero everywhere in 
the limit. In the region where (j) is bounded away from zero, the physical scalar curvature and the physical TT tensor 
remain bounded, while K 2 blows up. This is not compatible with solving the Hamiltonian constraint. Therefore <ft 
shrinks everywhere to zero, and the physical volume does likewise. 

We can do better. For ease of computation let us first make a conformal transformation so that the base metric 
has constant scalar curvature. The Lichnerowicz - York equation then becomes 



R - K% T K% T + ^K 2 = 0. 



(12) 




(13) 



8V 2 <50 - (8K% T Kj; T + \k 2 )5^ = \5K 2 . 



(14) 




(15) 



At the maximum of (j) we have V 2 (f> < 0. At this point we have 
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KZrKi? > Rv€ iax + ^K 2 ^ ax . (16) 

The K TT Ktt term is evaluated at the point where <j> reaches its maximum. We can replace this in inequality (16) 
by the maximum value over the whole manifold to give 

44 T H > R o€ iax + \K 2 ^ ax - (17) 

This means that as K 2 increases (f) max must decrease, and therefore the physical volume shrinks to zero. 

In the nonnegative Yamabe classes, where i?o > 0, one can just discard the i?o term in inequality (17) to get 

K^K^imax) > 2 -K 2 <f> 12 ax . (18) 

From this one gets 



v p = J ^dv < <tt ax v = y TT — V . (19) 

This estimate holds always in the positive and zero Yamabe classes; in particular, it shows how the volume goes to 
zero as K 2 becomes large. In the negative Yamabe class such an estimate will continue to hold, but only in the large 
K 2 limit. This can be seen by looking at the polynomial on the right hand side of (17) assuming Rq = —C 2 . 

The final question is what happens when K 2 shrinks to zero. While we know that the conformal factor increases 
and drags the volume with it, we have no guarantee that the volume becomes infinitely large. The Lichnerowicz - York 
equation, when K 2 equals zero, reduces to the original Lichnerowicz equation. Therefore we know that the Yamabe 
constant of the base metric must be the determining factor. If the base metric is in the positive Yamabe class, then 
the Lichnerowicz equation has a regular solution with finite physical volume, and a maximal slice emerges. This will 
form the limiting solution to the sequence of Lichnerowicz - York equations. The volume of the maximal solution is 
the largest volume we can specify. It, of course, depends both on the choice of base metric and on the specified TT 
tensor. 

If the Yamabe class of the base metric is nonpositive, then no regular solution to the Lichnerowicz equation exists. 
The conformal factor blows up as K 2 goes towards zero, and the physical volume becomes unboundedly large. 

Let us consider the case where the Yamabe constant is negative. This means that one can make a conformal 
transformation to a metric where the scalar curvature is constant and negative, i.e., R = —C = constant. One could 
pick the constant to equal —1, but this is not important. This base metric will have total volume Vq. As usual, one 
solves the Lichnerowicz - York equation for the conformal factor <fr, and calculates the physical volume V p . I will show 
that 

3(7 \*„ /3|Y|V 



The equality in Eq.(20) follows from the definition of the Yamabe constant, Eq.(9). This guarantees that the volume 
diverges as K 2 0. 

Start with a set of data, a metric with constant negative scalar curvature, a TT tensor, and K 2 . First make a 
conformal transformation with a constant conformal factor O = (3C/2K 2 ) 1 / 4 . This transforms the base manifold to 
one where the scalar curvature R = —2K 2 /3 and the volume V = (3C/2K 2 ) 3 / 2 Vn- The Lichnerowicz - York equation 
in this frame becomes 

8V 2 + -K 2 ^ + WTR^^ 7 - -K 2 ^ 5 = 0, (21) 
3 3 



where K tt Ktt — K TT Ktt(2K 2 /3C) 3 . From the min-max principle we know that at the minimum of <fi we have 



2- 

V 6 > 0. Therefore we have 



-K 2 ^ mm - min ) - KTT KTT ^ m > 0. (22) 

This guarantees that 4> min > 1 and thus <f> > 1. In turn, this means that the physical volume is greater than 
V= (3C/2K 2 ) 3 / 2 V . 
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We can do somewhat better than that by looking at the maximum of <j>. In this case we have 

\K\tnax 4> max ) K^K^^L < 0- (23) 

This can be rearranged to give 

^maAtmax ~ ^ ro aJ < \k^K^/K* = K TT K TT (2K 2 /'if 2_ . (24) 

This tells us that <f) max cannot be very large. As K 2 — » 0, the right hand side of Eq.(24) shrinks to zero and 
4> max is pushed down to 1. Thus the volume of the physical space converges to the volume of V = (3C/2K 2 ) 3 / 2 Vb = 
(3Y/2K 2 ) 3 / 2 . Another way of saying this is that in the limit, as K becomes small, the conformal factor becomes 
large, which renders the TT tensor unimportant. The physical geometry then more and more approximates a large- 
volume 3-manifold with constant negative scalar curvature R = —2K 2 /3. An alternativeway of expressing this is that 
V P K 3 w (3Y/2) 3 / 2 = constant. 

In summary: If the base metric is in the positive Yamabe class, and if we pick any physical volume which is less than 
or equal to the volume of the maximal solution, there exists a unique choice of K 2 , which, on solving the Lichnerowicz 
- York equation, generates a solution to the constraints with the chosen volume. If the base metric has nonpositive 
Yamabe constant, one can choose any volume between and oo; then there exists a unique choice of K 2 which delivers 
this volume, on solving the Lichnerowicz - York equation. 

III. OTHER CHOICES 

In addition to trading off between the physical volume and the value of K, there are several other specifications 
of data for the Lichnerowicz - York equation that work. Given a base metric, a TT tensor, and a constant K, one 
can find a complete family of TT tensors by simply multiplying the given tensor by a constant. Then, for each choice 
of the constant, one can solve the Lichnerowicz - York equation and find the physical volume. One can specify the 
volume and solve for the constant multiplying the TT tensor. 

A. Scaling the TT tensor 

We start off with a modified Lichnerowicz - York equation 

8V 2 - R<j> + a 2 K^ T K^ T <p- 7 - \k 2 ^ = 0, (25) 

and we let a 2 range from to oo. We now look for the relationship between the physical volume and the value of a 2 . 
Let us again perform a perturbation calculation, i.e., we start with a solution to the Hamiltonian constraint 

R-K^K« T + |fl* = 0, (26) 

and perturb Eq.(25) around it, holding the metric, the TT tensor, and K 2 fixed. We only perturb a 2 and consider 
the change of 4> away from 0=1. We get 

8V 2 <50 - R5(j) - 7K l ^ T K^ T S4, - ^-K 2 5$ = -Sa 2 K^ T K^ T . (27) 

By using Eq.(26) to eliminate R, we can simplify this equation to get 

8V 2 5<P - (8K% T K? T + rK 2 )60 = -Sa 2 K^> T K^ T . (28) 
o 

We can again use the min-max principle on this equation and if da 2 > we have 6<ft > 0. In this case, as distinct 
from the 'changing K 2 ' case, as we increase the K TT term the volume monotonically increases. This guarantees 
uniqueness. 

In order to investigate whether there is any restriction on the allowed range of the volume, we need to look at 
what happens at the two extremes, i.e., when a 2 — > oo and a 2 — > 0. Let us start with the case where a 2 — > oo. The 
Hamiltonian constraint, after solving the equation, is 
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R - a 2 r 12 Kj^K l S T + \k 2 = 0. (29) 

If (f> were to remain bounded on the support of K TT as a 2 — > oo, the central term in Eq.(29) blows up while the other 
two terms remain finite. This cannot happen; therefore <j> must blow up on the support of K TT , and the volume of 
the physical space becomes unboundedly large. This happens with any choice of base metric. 

The other extreme is when a 2 — > 0. In this case the TT term drops out of the Hamiltonian constraint, and we get 

R+\k 2 = Q. (30) 

This is a manifold with constant negative scalar curvature, and this is only achievable if the Yamabe constant of the 
base metric is negative. 

Therefore if the base metric is is in the negative Yamabe class, the conformal factor tends to a finite limit and 
generates a manifold of constant negative scalar curvature. This manifold of constant negative scalar curvature is 
the manifold of least volume in the sequence, and any physical volume we pick must be at least as large as this. We 
can find an expression for this minimum volume easily in terms of the Yamabe constant. On a manifold of constant 
curvature the Yamabe constant is achieved with a constant function. Therefore we have 



RV P 2K 



2 



V 3 



■V p * = Y. (31) 



From this we get 



This is the minimum volume we can specify if the Yamabe constant is negative. 

It is clear, from Eq.(32), that this mimimum volume shrinks to zero as \Y\ — ► 0. We can, however, do somewhat 
better. Let us assume that we are given a base metric in the positive Yamabe class. Let us make a conformal 
transformation so as to map the base metric to one with constant positive scalar curvature. Let us further arrange 
that this constant Rq = 2K 2 /3. This is not really necessary, but it makes the algebra somewhat easier. The L - Y 
equation now becomes 

8V 2 ^ - \k 2 4> + a 2 K TT K TT 4>- 7 - \k 2 ^ = 0. (33) 

At the maximum of 4> we have that V 2 </> < 0. This gives 

a 2 K TT K TT > 2 -K 2 {<f max + <^ 2 a J, (34) 

where, of course, the left hand side of Eq.(24) is evaluated at the point where <fr achieves its maximum. First, we see 
that 4> achieves its maximum on the support of K TT . Further, we can get an upper bound for <fi max in terms of the 
maximum of K TT Ktt from 

^(K TT K TT ) max > C- (35) 

This is not in any way sharp. A 'better' estimate can be obtained by turning inequality (34) into a cubic equation. 
This, however, offers no real improvement, because we are really interested in the case when a ~ * 0; since then 
(f>max — > and the <fi 12 term can be neglected wrt the </> 8 term. 
From (35) we get the following estimate for the physical volume 



V P = J 4?dv < 4L ax j dv = 4> & max V Q < 



(A- K TT ) r 



Vb, (36) 



where Vq is the volume of the base, constant curvature, metric. This shows that V p goes to zero as a goes to zero like 
c?l 2 . One could replace V in Eq.(36) by the Yamabe constant to give 



2K 2 



[{K TT K TT ) max Y . (37) 
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B. Scaling the base metric 



A third way of interpreting the Lichnerowicz - York equation is to choose a base metric g b , a TT tensor, a physical 
volume V p , and the trace of the extrinsic curvature K, but to assume the base metric is defined only up to a 
multiplicative constant g' b = fi A gb- The given TT tensor will still be TT wrt the adjusted metric. The value of (3 is a 
freely choosen parameter. 

The behaviour of this system exactly mirrors the second system analysed, where the TT tensor is scaled. This is 
because the following global scaling 

(g b , K) -> {(3 A g b , (3~ 2 Kj^ K) (38) 

leaves the physical metric, and also the physical volume, unchanged. Therefore there is a tradeoff between scaling on 
gb and a scaling of K TT . We have the following chain 

(g b , K? T , K) - ((3 4 g b , K% T , K) -> (g b , (3 2 K% T , K). (39) 

The first link is the nontrivial change we wish to analyse. The second link is an identity, which leaves V p unchanged. 
Therefore scaling the metric by (3 is exactly equivalent to scaling the TT tensor by 1 . 

This means that as (3 — > oo, the physical volume monotonically increases and becomes unboundedly large. As 
(3 — > 0, the behaviour depends on the Yamabe number of the base metric. If the Yamabe number is negative the 
volume shrinks to a finite volume, represented by a manifold of constant negative scalar curvature, while if the Yamabe 
number is nonnegative, the manifold shrinks to zero volume. The estimates given above in subsection A continue to 
hold, except that a has to be replaced by 1 . 
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